NASA TM X-67860 


NASA TECHNICAL 

MEMORANDUM 


N71- 27867 
NASA TM X-67860 



INTERFACIAL SHEAR FACTORS FOR LAMINAR 
FILM BOILING FROM A VERTICAL SURFACE 


by Kenneth J. Baumeister and Harold R, Jacobs 
Lewis Research Center 
Cleveland, Ohio 


TECHNICAL PAPER proposed for presentation at 
Thirteenth International Congress of Refrigeration 
Washington, D. C, , August 27 -September 3, 1971 



E-6381 


INTERFACIAL SHEAR FACTORS FOR LAMINAR FILM 
BOILING FROM A VERTICAL SERFAGE 

Kenneth J. Baiaineister 
NASA Lewis Research Center 
Cleveland^ Ohio 

and 

Haro Id R . Jac oh s 
University oT Utah 
Salt Lalce City;, Utah 

SLWiARY 

Forced flow film toiling heat transfer coefficients from a verrloal sux''f- 
ace have seen expressed in terms of the intex'’i acial shear or the so culled 
“slip" at the liquid vapor interface- The heat transfer coefficient is ex- 
pressed as a function of the Froude number and the shear fact or « The shear fac- 
tors are determined from an exact analysis v/hich accounts for the coupling he- 
tween the vapor and liquid momentiom at the interface. The correlating equa- 
tion may he useful to the designer interested in a simple correlating tecnnique 
valid over a wide range of Froude Numbers . 

INTRODUCTION 

Film boiling heat transfer plays an imporxant .role in energy transfer in 
cryogenic systems. In this paper film hoilin;- from a vertical surface will he 
studied. In particular^ w^e will study how the Froude n’omher and vapor ana 
liquid properties effect the interfacial shear or at the boiling inucr- 

face. A knowledge of the interfacial shear is important for high pressure cryo- 
gens .^specially in the vicinity of the thermodynamic critical point. 

The earliest heat transfer analysis of pool fimn boiling from a vertical 
surface was by Bromley (ref. l) . Later ^ Hsu and Westwater (ref. 2) ^./.tended 
Bromley's analysis for both laminar and turbulent vapor flo./. Bromley con- 
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sidered the equations of motion in the hot vapor film with the touna’;,ry c-.'a- 
ditions at the \7all equal to 

y = 0 u = 0 (l) 

Unfortunately, the houndary condition at the liquid vapor interface, c, sho'.m 
in Fiqure 1 case (a), was an unkno-wn. Bromley, however, considered two limit- 
ing possibilities on the velocity at the liquid vapor interface. The first was 
that (fig. 1(h)) 

y = 5g u = 0 (2) 

that is, the vapor does not drag any liquid along. The second possibility, 
considered was (fig. 1(c)) 

y = 6 m ^ = 0 (3) 

Sy 

that is^ the shear hetv/een the liquid and vapor is zero- Bromley derived a 
simple expression for the heat transfer coefficient for each of these tvc 
assessed boundary conditions. These expressions bracketed the unknown heat 
transfer coefficient. Bromley then suggested that a heat transfer coefficient 
which is the average of these two extremes be used as a correlating equation. 

Koh (refs. 3 and 4) extended Bromley's analysis by considering the momentum 
transfer between the liquid and the vapor. Koh's results are expressed nu- 
merically in terms of the important parameters such as Prandtl number and the 
density viscosity ratios. From his analysis^ Koh showed that for 






< 0.01 
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■boundary condition (eq. (2)) is a good choice. Equation (4) generally ho ..cis 
for low pressure boiling. A notable exception to equation (4) is film boiling 
near the thermodynmic critical poinu. 

Later forced flow past a vertical surface in film boiling was studied 
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(figs. 1(d) and (e)). Rankin considered the equations of raouion in the va,-- .r 
and av the liquid vapor interface, he assumed 

y = 6 u = u„ (3) 

that is, the interface moves at the free stream velocity. For high liquid 
velocities this is a reasonable assumption if (see ref. 6). For a 

low liquid velocity or near the critical point, however, this is only an 
approximation as seen from the flow profile shown in Figure 1(d) . 

Recently Jacobs and Boehm (ref. 6) extended the theory for the case of 
film boiling from vertical surfaces with flow. They coupled the governing 
equations of motion and energy in both the vapor film and surrounding liquid. 

In this case, the numerical solutions of the heat transfer coefficients are 
correlated as a function of Froude number as well as other system parameters 
included in Koh's (refs. 3 and 4) pool boiling analysis. The nmerical results 
in reference 6, however, are limited to a single curve at one density- 
viscosity ratio. 

The purpose of the present paper is to present a relatively simple heat 
transfer correlation for flow film boiling past a vertical surface. First, a 
shear parameter K is defined; then, the equations of motion and energy in the 
hot vapor are solved with the factor K treated as -an unknown parameter. 

Finally, by a comparison with the exact numerical results in reference 6 the 
actual K is specified. The simple heat transfer correlation which results 
may be of practical use to the design engineer. 

LIST OF SYMBOLS 

Op specific heat at constant vapor pressure 

Fr Froude Number u^/gx 

Or Grashof Number py(pj^ - py)gx^/n2 

g coefficient of gravity 
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Sc 

h 

h* 

K' 

Nu 

P 

Pr 

S 

Tw 

Ts 

AT 


gravitational constant 

heat transfer coefficient 

dimensionless h defined "by equation (12) 

shear factor^ equation (lO) 

thermal conductivity of vapor 

Nusselt number hx/ky 

pressure 

Prandtl number jiCp/ky 
defined by Cpifli/APr 
wall temperature 
saturation temperature 

dw - Ts) 


u 

^max, pool 
W 



l-^V 

a 


velocity parallel to surface 

free stream liquid velocity 

see equation (ll) 

mass flow rate per unit length 

distance along plate from leading edge 

distance normal to surface 

vapor gap thickness 

vapor gap thickness for complete slip (K = O) 

vapor gap thickness for no slip (K = l) 

latent heat of vaporization 

viscosity of liquid 

viscosity of vapor 

surface tension 


Theory 

The following analysis for laminar film boiling along a vertical surface 
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assumes a smooth laminar film^ constant pressure Gradient^ constant vapor prop- 
erties which are evaluated at the film temperature^ and negligihle fluid inertia. 
In the energy analysis) heat is assumed to he transferred only by molecular con- 
duction. The liquid is assiamed to be at its saturation temperature. Conse- 
quently, the governing momentum and energy equations for the vapor film shown 
in Figure 1 are as follows: 



5^u ^c 
By2 PV <5x 


a = 

W = 



h AT 


dy 


(5) 

( 6 ) 
(7) 


Equation (s) is the momentum equation; equation (6) is the energy equation; 
and equation (?) accounts for the increased mass vapor flow due to vauporization. 

, The boundary conditions on the energy equation are simply 

y=*0 T = T^;y=6 T = Tg^^. (8) 

The boundary condition on the momentum equation at the wall is given by equa- 
tion (l). The boundary condition at the liquid vapor interface, however, 
needs special consideration, and it is here that we introduce the slip ratio. 
Shear Factor K 

For pool boiling^ the no shear boundary condition^ equation (3)^ leads to 
a predicted vapor velocity which is the maximiom possible obtainable at the 
liquid vapor inoerface^ "^max^pool* However^ in ilow boiling at high Froude 
numbers the free stream velocity u^ controls the interfacial velocity. Con- 
sequently^ because of interfacial shear^ the actual interfaciaC. vapor velocity 
is within the range 


0 < u(5) < u^ + Uj^a^^pooi 


(9) 
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For high Fr momhers, dominates while at low Fr numhers or pool ho-^ling 
^ax nool dominates. From these considerations, the choice of a houndary 
condition at the liquid vapor interface would he 


y - 6 u - K(uoo + 'Umax, pool) (lO) 
where K is a shear factor which depends on the magnitude of Uoo, fluid prop- 
erties, and position x. It can he readily shown (ref. l) that 


( 6 ) _ ~ 
^ax,pool “ 


■kY^fl?Xp,v 


hy L^Py^PL ” Py)sJ 


1/2 


( 11 ) 


Solution for the Heat Transfer Coefficient 

The solution of equations (5), (6), and (7) subject to boundary conditions 
(8), (l), and (lO) yields the following expression for the film boiling heat 
transfer coefficient written in terms of Froude number and the factor K 
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Because of paper length restrictions, the derivation is not given here. How- 
ever, the integrations are simple and the integration steps follow in a simi- 
lar manner to that given in reference 5 for the no slip case. 

In the limiting case of pool boiling (u„ = 0 or Fr = 0), for maximum 
shear (K = O) h* = 0. 5, while for no shear (K = l) h* = 0. 707. These results 
agree with the earlier work of Bromley. 

Solving equation (12) for the factor K yields 
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DISCUSSION 
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In reference 6. a single curve is presented for a dimensionless rx-at 'crxis-- 
fer coefficient as a function of Fr number for S = GxlO”*"^^ 

I^l/uy = 25. A plot of the factor K for this curve is shown in Figure 2 
as function of Froude niomber. As seen in Figure 2^ for Froude numbers less 
than 0.001 the K parameter approaches a constant value of 0.59;, while for ' 
Froude number greater than 100^ K reaches a value near 0.75. 

Near a Fr number of 0.01 a minimum in K is seen. This dip could be due to . 
the manner in which K is defined in the low velocity region;, where u^ and 
%ax^pool same order. 

The asymptotic value of K at large Froude numbers for this case was 0. 75;, 
because of the small value of S assumed in reference 6. In general;, however^ 

K will approach an asynptotic value of 1 for large Froude numbers whenever the 
inequality given hy equation (l) hold true^ as shown in Figure 2. 

At the present time more exact numerical solutions are required to map K 
over a large range of density and viscosity ratios especially near the thermo*- 
dynamic critical point where the density- viscosity ratio approaches one. Next;, 
a proper set of dimensional groups should he derived such that all the curves 
could he reduced to a single curve. 
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